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We show that both fermionic and bosonic uniform d-dimensional lattices can be reduced to a set 
of independent one-dimensional chains. This reduction leads to the expression for ballistic energy 
fluxes in uniform fermionic and bosonic lattices. By the use of the Jordan-Wigner transformation 
we can extend our analysis to spin lattices, proving the coexistence of both ballistic and non-ballistic 
subspaces in any dimension and for any system size. We then relate the nature of transport to the 
number of excitations in the homogeneous spin lattice, indicating that a single excitation always 
propagates ballistically and that the non-ballistic behavior of uniform spin lattices is a consequence 
of the interaction between different excitations. 
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INTRODUCTION 

Due to the rapid advance in nano-technologies the study of non-equilibrium transport phenomena in quantum 
systems, including charge and heat transport, has become a major field of study. Furthermore, transport properties 
in spin and harmonic oscillators lattices are important because they model many realistic physical systems. Spin 
lattices and ladders can be experimentally realized by different methods m- Harmonic oscillators lattices describe 
a plethora of different systems including optical cavities arrays [4], trapped ions [5], and phonons in cubic crystal 
lattices iniiz]- Finally, fermionic and bosonic atoms transport can be studied in optical lattices showing both ballistic 
and diffusive regimes pin] 

The role of dimensionality is crucial in both quantum and classical systems |10j . One-dimensional systems composed 
by spins mm and harmonic oscillators m have been previously analysed by the use of Markovian master equations. 
These analyses show that these systems are ballistic when the transport takes place coherently and become diffusive 
if a dephasing channel is locally coupled. This is a typical behavior of one-dimensional homogeneous systems in the 
presence of noise, while disordered systems possess additional features M- The same behavior is found in harmonic 
lattices of arbitrary dimension |18j . However, in multidimensional spin lattices the results are different, and an 
analytical solution is lacking. The simplest two dimension topology, a ladder, has been analysed by Znidaric m, 
showing that quantum transport is in general anomalous. In fact, this system has a set of ballistic subspaces, making 
it possible to create ballistic transport by the design of non-local baths operators and the initial state preparation. It 
has also been demonstrated experimentally with ultracold atom in optical lattices that spin systems of one and two 
dimension behave very differentlyj^. 

Another canonical problem regarding multidimensional transport is quantum walk (QW). It is known that a QW 
in a one dimensional system propagates faster than classical random walks m- In multidimensional systems, it has 
been proved by Kempe m that a discrete QW in a hypercube hits from corner to corner exponentially faster than 
a classical random walk. This result has been extended to distorted hypercubes [18] and to hypercubes embedded 
in more complex graphs m- Finally, in hypercubic lattices some searching algorithms based on QW have been 
developed showing a speed-up against their classical counterparts noiin]- These results suggest that qubit lattices 
with only one excitation behave in a superdiffusive way. 

Regarding the transport in quantum lattices there are still some open questions: 

• What is the role of the fermionic and bosonic nature of the system in the transport properties? 

• Why spins and harmonic oscillator systems behave differently in lattices with dimension d > 1? 

• What is the origin of the ballistic and non-ballistic subspaces found in spin ladders? 

• How does the number of excitations in a spin lattice affect transport properties? 

In this paper we study the general case of quantum transport in nonequilibrium lattices composed of fermionic or 
bosonic sites. We show that these lattices with an arbitrary dimension can be decomposed into independent one¬ 
dimensional normal modes. This confirms that both fermionic and bosonic lattices are ballistic in any dimension and 
allows us to derive explicit expressions for energy fluxes of uniform fermionic and bosonic lattices. We then use this 
method to analyse spin d-dimensional lattices showing that ballistic spaces exist in any dimension and size of the 
system. Finally, we discuss the implications of this result in hypercubic lattice quantum walks and discuss the role of 
the number of excitations in transport properties. 


LATTICE MODEL 

The systems we study are fermionic and bosonic d-dimension lattices. In order to analyse the features of the 
transport we couple the system to incoherent thermal baths that drive it out of equilibrium. Due to the effects of the 
baths, the system evolves to a time-independent steady state. If the temperatures of the baths are different, there is 
a finite energy current flowing through the system. The dependence of the current on the system size indicates the 
nature of the transport. We consider transport to be ballistic if the current is independent of the system size mis]. 

In figure [T] a d = 3 lattice with local thermal baths at the end of dimension 3 is displayed. Each site of the 
lattice is defined by d indices li = 1,..., Lj, where Li is the size of the system in dimension i. We define the site 
vector 1 = {li,... ,ld) and the size vector L = (Li,..., Ld). The lattice is connected to thermal baths at the ends 
of dimension d, making this the direction where the energy transfer takes place. Periodic boundary conditions are 
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FIG. 1. Sketch of a d = 3 lattice with L\ — 3, L 2 = 2, and L 3 = 3. Local thermal baths apply at the terminal sites of dimension 

3. 

applied to all the remaining dimensions. We also define the reduced vectors to represent all the dimensions but the 
last one, 1^ = {li,..., Id-i) and = (Li,..., L^-i). 

The full Hamiltonian of our system can be decomposed into the free Hamiltonian, which depends only on the 
occupation number of each site, and the hopping Hamiltonian in each dimension, which is responsible for transport. 
We can write the total Hamiltonian as H = iJfree + h=l throughout the paper) 


^iree — ^ ^ ^ 

1 

Hi = (“L- ,ii,- .h+1.- ,id + ’ 

1 


( 1 ) 


were w is the frequency of each site, and gi represents the coupling strength in dimension i. The operators a and are 
ladder operators that can be either fermionic (/ ,/^) or bosonic (6 ,b^), obeying anticommutation and commutation 
relations respectively. In the fermionic lattice, as there are no internal degrees of freedom each site can contain at 
most one fermion (//// = 0). On the contrary, the bosonic lattice can contain an arbitrary number of bosons per 
lattice site. 

The system is driven out of equilibrium by bosonic thermal baths locally coupled to the end of the system in 
dimension d. The overall dynamics of the system is described by a Markovian master equation [22j . 


p = -i[H, p\+ Cip + Cl^P 

Here, each of the reservoirs is modelled by a Lindblad super-operator 


( 2 ) 


+ X! ^'^/Ld (^1/id + 1) (%,l2,--;l/LdP “LL....,l/id “ 2 {“h.L,....l 

Ir- ^ 






( 3 ) 


where £1 acts at the beginning of dimension d and models the bath attached at the end of this dimension. The 
first term in Ci accounts for emission into the reservoir, and the second term for absorption, T is the interaction 
rate, and Ui is the mean excitation number at the resonance frequency of the reservoir. For instance, in the case 
of a bosonic bath at temperature T^, rii = l/[exp(a;/Ti) — 1] (with Boltzmann’s constant ks = 1). We perform the 
analysis using Bose-Einstein statistics, but the results are independent of the nature of the baths. 
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The validity of this master equation in nonequilibrium systems has been extensively analysed. In Ref. |23] a system 
of two interacting harmonic oscillators coupled to thermal baths at different temperatures is studied both numerically 
and analytically, concluding that the master equation ([^ is valid in the limit of small intercoupling. Furthermore, in 
Ref. [13] it is proved that in the case of bosonic lattices with two local baths at the same temperature the system 
thermalises in the limit of small intercoupling. 

The energy transfer through the lattice can be quantified by the time derivative of the energy expectation value in 
the system 


E=j^{H)=Tt{Hp). (4) 

This expression should be zero at the steady state. Using the master equation (|^ we can decompose Q into 

E = Tr (HLip) + Tr {HC 2 P) := Ji + J 2 = 0, (5) 

where Ji and J 2 are equal in magnitude but with opposite signs. Each of these terms refers to the mean energy 
interchanged with each thermal bath per unit of time. This allows us to define the energy flux through the system as 
J = I Ji I = IJ 2 1. In regular systems the behaviour of the energy flux is proportional to the behaviour of the spin or 
excitation flux. In more complicated systems like networks the energy and excitations fluxes can behave differently 

m- 


NORMAL-MODE DECOMPOSITION 


To calculate the energy transfer through an arbitrary lattice, we define normal modes associated with all the 
dimensions except for d, q^°‘ = where a = 1,... ,d — 1 and Ua = I,..., Lq [331 US] ■ Therefore, we define the 

mode vector as q = and there are N = nf=i normal modes. We transform the ladder operators 

to the normal mode basis via the following transformation 




a 


t 


1 

^/N 


q 


d-l 

27rz^ 

i=l 



1 

^/N 


q 


d-l 

2=1 



( 6 ) 


where and a\^ are ladder operators in the normal mode basis that fulfils the same commutation/anticommutation 
relations as the site ladder operators. Applying this transformation to the Hamiltonian Q decomposes it into 
independent terms for each mode, H = X]q^(*l) with 


/ d—l \ Ld Ld—i 

H{q) = ( w + 2 ^ gfc cos (g”“) ) JI “L + 9d Y (“L “'d+i 

V fc=i / /d=i id=i 

Also the Lindblad super-operators <§ decompose into independent terms for each normal mode 
Cl/Ld (q) P = r l/Ld^l/Ld ((^l/Ld P°-\/Ld (q) “ \ {“I/Ld “l/Ld > p}) 

+ Ti/id {ni/Ld + 1) (al/Ld (q) P“i/Ld (q) {^I/Ld (q) “I/Ld (q) ’P}) ■ (8) 

The master equation ([^ unravels into a set of N independent equations, each of them corresponding to one of the 
modes, which behaves as bosonic or fermionic one-dimensional system. Each of these equations is equivalent to the 
equation of a one dimension system with the on-site energies shifted by 2 9k cos (g2“)- 

Hence, the problem of calculating the energy flux through a d-dimensional lattice can be reduced to calculating the 
energy current in a one-dimensional chain. The heat transfer through one-dimensional chains composed by harmonic 
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oscillator has already been solved in Ref. |13j . The one-dimensional fermionic chain is solved in the Appendix. The 
analytical expressions for the energy flux for each mode is 


Jf (q) 
(q) 


d-l 


w + 2^gfcCOS (g2“) 

V / 

/d-l \ 

w + 2^gfeCos (g2“) 


/c=l 


4gd7i7Lj(si - 
(71 + 7L J (4gd + 7 i7l J 

^9d^i^LAni-nL^) 

(ri + rij(4g2 + rir,.j’ 


(9) 


where we have introduced the notation 7 ^ = Ti{2ni 1) and Si = ni/(2ni -f 1) in the fermionic case. The total energy 
flux can be calculated by summing up the flux of each channel, Jtotai = N Jij_ where Jid is the energy flux of a uniform 
one-dimension system with on-site energies uj. 

The lattice size in the energy transfer dimension plays no role in the amount of energy transfer, as it does not appear 
in This result shows that both fermionic and bosonic transport in the absence of noise or disorder are ballistic 
with arbitrary size and in any dimension. As our model coincides with the Hubbard model with no interactions 
(U=0), this result explains why in the non-interacting regime there is ballistic transport in optical lattices both with 
fermionic and bosonic atoms liH]. 

Even though both fermionic and bosonic transport are ballistic, their dependences on parameters such as the 
temperatures are different. In the limit of zero temperature the two converge as limr.-^o 7 i = T^ and limr^-^o Si = rii. 
But the situation at finite temperature is very different. The energy currents of one dimensional fermionic and bosonic 
systems are displayed in figure as functions of the temperature of the hot bath (the cold bath is kept constant at 
a very low temperature). For small values of the hot bath temperature the currents are very similar in magnitude. 
At higher temperatures the bosonic system has a higher energy current than the fermionic one. This is a consequence 
of the Pauli exclusion principle that limits the heat capacity of the fermionic system. In the limit of very high 
temperature of the hot bath the current of the fermionic lattice decreases as the temperature increases, similar to 
spin systems [ 12 ]. 



1 2 3 4 5 


T 2 


FIG. 2. Energy current (logarithmic scale) for a fermionic and bosonic system as a function of the temperature of the hot 
bath. Fi = r 2 = 0.01, g = 0.01, uo = 10, Ti = 0.001, d = 1. 

Finally, we need to point out that there are also non-ballistic systems that can be decomposed into non-interacting 
one-dimension channels by the application of normal mode transformation ([^. The mode decomposition can be 
applied to a broader set of Hamiltonians than defined in Q. If we relax the condition of homogeneous frequencies 
and couplings to the condition of having them equal in all the dimensions except for the direction of transport, then 
the inhomogeneous system can be decomposed into one-dimension disordered chains. These disordered one-dimension 
systems have been broadly studied and are known to display numerical evidence of non-ballistic behaviour [21 [55] . 








6 


SPIN LATTICES 


In one dimensional spin systems driven out of equilibrium by thermal baths (XY model) the transport is similar to a 
homogeneous fermionic chain. In the spin case the current is given by an expression similar to (|^ |12j . This is a natural 
result as one dimensional spins systems can be mapped locally to fermions by the Jordan-Wigner transformation m- 
Spin chains with different Hamiltonians have been broadly studied nmnnH], concluding that in most of the cases 
the transport in a homogeneous and noise-free spin chain is ballistic and it becomes diffusive in the presence of static 
or dynamical noise. 

On the other hand, the similarity between spin and fermionic systems breaks down if the dimensionality is larger 
than one. As we have proven, fermionic lattices are ballistic in any dimension as long as the couplings are homogeneous 
in each dimension. But it has been demonstrated that this behaviour does not hold for spin lattices. Even the simplest 
two-dimension spin lattice, a homogeneous ladder, presents both ballistic and non-ballistic invariant subspaces. The 
relative size of the ballistic subspaces goes to zero as the size of the ladder goes to infinity [15]. This fact highlights 
the difference between spins and fermions as well as the importance of the dimensionality in spin transport. 

A general spin lattice of dimension d has a Hamiltonian similar to the fermionic lattice Q. It can be decomposed 

in the form //■!»"■ = i/£“ + H'""" with 


1 


( 10 ) 


where ja are the Pauli raising/lowering operators. 

Due to the similarity of Hamiltonians Q and (10) it is natural to think that the mode transformation Q should 
also decompose a d-dimension spin lattice into independent one-dimension spin chains, but this is not the case. Both 
bosons and fermions have uniform commutation and anticommutation relations respectively, but spin systems have 
not 0 Therefore, by applying the transformation Q to a d-dimension spin lattice, one-dimension mode operators 
can be defined. But these operators have different commutation rules compared to the original operators in the site 
basis. Because of that, this transformation does not map a d-dimension spin lattice into independent one-dimension 
spin chains. 

Nevertheless, spin systems can be transformed to fermionic systems by the Jordan-Wigner transformation m- 
This transformation is dehned by 


fi = n n 

j<i j<i 

where the index i refers to an arbitrary ordering of the spins. Owing to this arbitrariness there is not a unique 
representation of the Jordan-Wigner transformation. In one-dimension systems with next-neighbours couplings this 
transformation can be performed in such a way that a homogeneous spin chain is mapped into a homogeneous fermionic 
chain. This explains why these two systems have the same behaviour. On the other hand, in systems with more than 
one dimension this is not true. Due to the non-local character of the Jordan-Wigner transformation, when it is applied 
to a multidimensional spin Hamiltonian the corresponding fermionic Hamiltonian includes also non-local terms. 


^ Note that if we have a spin system with sites labelled by the index i we can define the dynamics by the Pauli operators i 'rf |. 

We can then define a set of pseudo fermionic operators by stating crd = /I, af = f^ , and erf = 2fjf^ — I. These operators fulfil 
the anticommutation relation J- = 1, but operators from different sites commute |^/l,/j j = 0 (Vi yf j). There is no local 

transformation that transforms the Pauli operators into fermionic operators. 
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FIG. 3. (Left) Spin ladder coupled to two thermal baths at different temperatures. (Right) Labelling of the spin ladder. 

By the use of this transformation, we analyse the simplest spin lattice beyond one-dimension, a ladder. We order 
the sites by defining a new index I = li + 2(^2 — 1), where = 1,2 and I 2 = I, ■ ■ ■, L, and L is the size of the system 
in the direction of the energy transfer (see figure]^. The Hamiltonian of a uniform spin ladder in terms of the Pauli 
matrices is 


2L 


2L-2 


^iTddL = + H.c.) +92^2 ^1+2 + ■ 

1^1 1^1 1^1 

After applying the Jordan-Wigner transformation 0 we obtain the fermionic Hamiltonian 


( 12 ) 


2 L 


-^ladder — ^ ^ ^ // // 


9i 


1^1 


L 

E 

1^1 


(/2i/2Z-l 


2L-2 


-bH.c. 


92 




+ H.C. W+1 


1^1 


(13) 


where W = fj. This expression is similar to the fermionic ladder given by U , but it includes a new term 

that is responsible for the non-ballistic behaviour already found in this system. On the other hand, as stated in Ref. 
m this system has a plethora of invariant subspaces. By direct analysis of Hamiltonian 0 the transport properties 
of each subspace can be inferred. If any state within a subspace, \k), fulfils that Ni+i\k) = c\k) with a constant value 
of c for those values of I that satisfy ///;_|_ 2 l^) 7 ^ 0 or f 22 fi 1 ^) ^ 0 ; system reduces to a uniform fermion ladder 
and the transport is therefore ballistic. 


The role of the number of excitations in the non-ballistic character of the system is clear from (13). If the system 


contains only one excitation it is equivalent to a uniform fermion ladder, as in this case, and the transport is ballistic. 
This result can be extended straightforwardly to a more general case of a d-dimension lattice with arbitrary size, 
proving that there can be a ballistic spreading of the excitation in all directions and extending previous results about 
discrete quantum walks in high dimensions |18H20] to the continuous domain. 


Examples 

Invariant subspace 

In Ref. [15| Znidaric designed an open spin ladder with ballistic transport. This design is based on the use of 
non-local Lindblad superoperators that keep the system in a ballistic subspace. Together with a proper initial state 
he proved by numerical simulation that the flux is independent of the system size. The analysis of this system is 
simplified by using a Bell-type basis in the ladder rungs. This basis is defined by the vectors 

\S)i = ^ {crt +<^t+i) |vac),,,+i, 

1 ^)* = |vac),,,+i, 

1^)2 — 

\I)^ = cr+cryi|vac)*,,+i. 


( 14 ) 




with i being an odd number and |vac)i^i+i being the vacuum state of sites i and i + 1 . 

The initial state of the system is composed by alternating entangled states in the form \STST ...) or \TSTS ...), 
and the Lindblad superoperators act non-locally 


Li = \ITS ...ST){STS ...ST\, 

L2 = \IST ...TS){T ST ...TS\, 

L 3 = |S'TS'...S'r)(S'TS'...S'/|, 

U = \ST ...TS){TSTS (15) 

The effect of the baths in this case is given by the superoperator 


4 


(16) 


We have not specified temperature or coupling strength for simplicity, as the ballistic character of the system is 


independent of these parameters. By applying the Jordan-Wigner transformation (111 the Bell basis elements can be 
written as a function of fermionic operators. 


\S)^ = + /Ifi) |vac)i,*+i, 

\T)^ = - fhi) |vac)i,*+i, 

|0)j = |vac)i,i+i, 

l-f)* =////+i|vac),,i+i. (17) 

The ballistic character of this subspace can be proved by direct inspection. All the states in the form \ST)i (or 
\TSh, i+i) are eigenvectors of the Hamiltonian and they are therefore invariant. Furthermore, for all the states \k) in 
the form (or ) with ///;_,_ 2 |fc) ^ 0 or \k) ^ 0 one obtains Ni+i\SI) = \SI) 

due to the entangled state of the background. Consequently, this subspace is ballistic, as it can be mapped into a 
homogeneous fermion system. 


Double exciton subspace 

In the absence of hot and cold baths, if a system is prepared in an initial state that belongs to a ballistic subspace 
it will propagate ballistically without leaving the subspace. To study this behavour we have simulated the dynamics 
of a closed ladder. This allows us to analyze different initial conditions in a way that is impossible with a system 
coupled to baths. First, we study the dynamics of a spin ladder with L = 130 and gi = 92 = 9 with an initial state 
given by 


IV'(O ))0 = + e^*(ytcrt+2) \^ac), (18) 

which is a superposition of two excitations on the top and bottom legs at the middle of the ladde as it is displayed 
in figure 4 (left)r. To analyze the transport properties we have calculated the second derivative of the mean square 

displacement of the excitations in the direction of the energy transfer, given by C = ^ . This is a common measure 

of transport properties |29]. For ballistic transport we have a constant value C = where n is the total number of 

excitations in the system. The results are displayed in figure ([^ for different values of </>. Only for (f> = tt the system 
is in an invariant subspace corresponding to a uniform fermion ladder and the transport is ballistic. This prediction 
is confirmed in figure right by the constant value of C = 8 . The opposite behaviour is found with (j) = 0, and it is 
characterized by strong oscillations of C. 
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FIG. 4. Left: ske tch of the initial state (181. Right: C for a ladder as a function of time for different values of (j). The initial 
state is given by (18l. 


Four exciton subspace 

We provide another example inspired by Ref. m- We first define the Bell-type state 

l^0)i = ^ + e*'^cri+i) \vac)i^i+i. (19) 

It is clear that |^ = 0) = l^) and |^ = tt) = \T). The initial state is defined as 

\ipm^ = \O...OSA^SA^O...O), (20) 

where only the central sites are excited to alternating S and A^ states. With (f) = t: one restores the ”S — T 
background” introduced in Ref. where all other sites are in the ground state. The results with four excitations 
are shown in figure ([^. Similar to the previous example, it is proven that with (p = n the system stays in the ballistic 
subspace indefinitely. This can also be generalised to states with arbitrarily higher or odd numbers of excitations. 
For example, the state |'!/'(0)) = \0... OSTSO .. .O) containing three excitations belongs to a ballistic subspace. 



FIG. 5. G(t) for a ladder as a function of time. The initial state is given by (201 with four excitations in the center of the 
ladder. We follow the same convention used in hgure 
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All the ballistic subspaces of the spin ladder described in Ref. |15j correspond to uniform fermion ladders. This 
result can be easily extended to general lattices, proving the existence of invariant ballistic subspaces for a spin lattice 
of an arbitrary size and dimension. For an open system the transport properties depend on the initial state and on 
the explicit form of the Lindblad superoperators that model the interaction with the baths, so that it is possible to 
design systems with multiple steady states with different currents [3D]. This allows the design of a spin lattice with 
ballistic and non-ballistic steady states and the control of the system current by regulating the symmetries of the 
system. 


CONCLUSIONS 

We have presented in this paper a normal mode decomposition scheme that reduces d-dimensional bosonic and 
fermionic uniform lattices of arbitrary dimension into independent one dimension chains. These lattices are relevant 
for the study of many realistic systems such as optical lattices and cubic crystal lattices. By the use of the normal 
mode technique, the analytical expression for the energy transfer in an arbitrary uniform lattice is derived. For 
any dimension and size of the lattice, and with local heat baths acting at the ends of one of the dimensions, the 
system behaves ballistically and the energy transfer flux is independent of the system size. This result holds for both 
fermionic and bosonic systems. However, the dependence of the energy transfer on the bath temperature is different 
for fermionic and bosonic lattices, especially at high temperatures. 

This method can also be extended to spin lattices, which has been previously proven as non-ballistic. The analysis 
is carried out by applying the Jordan-Wigner transformation, that transforms a spin system into a fermionic one. 
The application of mode-decomposition to spin lattices proves that in general a uniform spin lattice corresponds to a 
non-uniform fermionic lattice. Spin lattices have many invariant subspaces due to the high symmetry on the system 
topology. Some of these subspaces are equivalent to fermionic uniform lattices and are consequently ballistic. The 
remaining subspaces are equivalent to non-uniform fermionic systems. This result generalises previous studies of spin 
ladders to systems with arbitrary size and dimension. 

Finally, by applying the mode-decomposition technique we demonstrate that a single excitation in a d-dimension 
spin lattice always propagates ballistically. This observation is relevant for the development of algorithms that use 
continuous quantum walks in a hypercubic lattice. 
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Appendix A. Calculation of the energy transfer in a fermionic chain 

The calculation of the energy transfer in a one-dimensional fermionic chain can be performed by applying the same 
method used in Ref. m- The one dimensional system Hamiltonian reads 

L L-l 

were ////^ are the creation/anhilitation fermionic operators, L is the number of sites of the system and h = 1. 

The system is driven out of equilibrium by bosonic heat baths connected to the terminal qubits and describe by 
the Lindblad terms 


^\/lP 


= hlLisilL — 1 ) (fi/LPfl/L ~ 2 



(A.22) 


The full dynamics of the system is given by the equation 


p = -i [H, p] + Cip -I- ClP- 


(A.23) 


The energy flux is given by 


J = Tv{HCip) 


^ 71 (si - {flh)) - ^ ((/ 1 V 2 ) - (a/I)) 


(A.24) 


with = Ti{2ni -I- 1) and Si = -f 1). The heat current at the steady state can be calculated only from the 

excited state population of the first site and the imaginary part of the coherence between this site and the second one. 


Due to the structure of the master equation (A.23) the next-neighbouring sites coherences are purely real, and the 
energy flux depends only in the first site population (see Ref. [H]). Therefore, the problem of calculating the energy 
flux is reduced to the problem of calculating the population of the first site at the steady state. 

At the steady state the time-derivative of each expectation value is equal to zero. By calculating the factors 
S (/fc/fc) ~ fc = 1, L we obtain the pair of equations 


71 (si - (/iVi)) = -9 ((/ 1 V 2 ) + (/i/l)) 

7L (sL - (/1 /l)) = 9 [{iI-Jl) + (/l-i/I)) > 

(A.25) 

for the terminal sites. Whereas for the inner sites the result is 

(fl-ifk) + (fk-ifl) = (/fcVfc+i) + (fkfUi) • (A.26) 

The last equation is obtained by summing up the coherences, ^ X)fc=i (^flfk+i^ ~ 

y (4/2) + Y (fl-ih) = -^9 ((/iVi) + {flh)) (A.27) 
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By solving these five equations, we can calculate the first site population 
flux equation (|A.24|) we obtain 


flfi 


J 


4 g^ 7 i72(51 - 52) 
(71+72) ( 4 g^ + 7172) 


and by substituting in the energy 

(A.28) 




